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Radiative annihilation of a soliton and an antisoliton in the coupled sine-Gordon 
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In the sine-Gordon equation solitons and antisolitons in the absence of perturbations do not 
annihilate. Here 1 present numerical analysis of soliton-antisoliton collisions in the coupled sine- 
Gordon equation. It is shown that in such a system soliton-antisoliton pairs (breathers) do annihilate 
even in the absence of perturbations. The annihilation occurs via a logarithmic-in-time decay of 
a breather caused by emission of plasma waves in every period of breather oscillations. This also 
leads to a significant coupling between breathers and propagating waves, which may lead to self- 
oscillations at the geometrical resonance conditions in a dc-driven system. The phenomenon may 
be useful for achieving superradiant emission from coupled oscillators. 

PACS numbers: 05.45.Yv 74.78. Fk 42.65.Wi 05.45.Xt 



INTRODUCTION 

Analysis of soliton dynamics in the sine-Gordon (SG) 
formalism is important in many research areas [1], includ- 
ing non-linear optics, condensed matter physics, atomic 
[2] and particle physics [3". Solitons are elementary par- 
ticles of the sine-Gordon equation, in a sense that they 
are quantized and do not spontaneously decay. Soliton- 
antisoliton collision is a nontrivial example of interaction 
of strongly non-linear waves. It may lead either to pas- 
sage of the two waves or to formation of a bound pair 
- the breather [4]. Within the pure SG equation soliton 
and antisoliton do not annihilate because the collision 
is elastic and the annihilation is prohibited by energy 
conservation. However, addition of various perturbation 
terms to the SG equation does allow particle-antiparticlc 
annihilation via breather decay [H |5]. This may hap- 
pen both via intrinsic viscous damping and via radiative 
losses from the breather [1 . 

In recent years properties of solitons in the coupled 
sine-Gordon equation (CSGE) are being actively studied. 
The CSGE describes complex behavior of interacting sys- 
tems, such as atoms in a periodic potential [5], magnetic 
multilayers [51, stacked Josephson junctions [7HTU] and 
layered superconductors ^llril3i . Coupling of N systems 
leads to a variety of unusual effects. First of all, it leads to 
appearance of N eigenmodes with different symmetries, 
length scales and velocities [51 [Uj . Even though the ex- 
act soliton solution in this case is not known, numerical 
and approximate analytic results demonstrated that the 
soliton becomes composed of different eigenmodes [9l [10] 
and the shape of such a composite soliton may become 
very unusual in the dynamic case. Next, unlike the sine- 
Gordon equation, the coupled sine-Gordon equation is 
not Lorentz-invariant |13] . Therefore superluminal soli- 
ton motion (faster than the slowest eigenmode velocity) 
is possible |9l [131 IE] ■ It is accompanied by Cherenkov- 
type radiation, due to decomposition of soliton compo- 
nents with eigenmode velocities slower than the speed of 



the soliton into plasma waves travelling along with the 
soliton |9l[ini[l5]. 

In this work I present numerical analysis of soliton- 
antisoliton collisions within the coupled sine-Gordon 
equation with focus on Josephson vortex (fluxon) dynam- 
ics in magnetically coupled stacked Josephson junctions. 
Both direct (fluxon and antifluxon in the same junction) 
and indirect (fluxon and antifluxon in different junctions) 
collisions are considered. It is demonstrated that soliton- 
antisoliton pair in the CSGE can annihilate even in the 
absence of viscous damping or other perturbations. An- 
nihilation occurs via emission of plasma waves from an 
oscillating breather, to some extent resembling annihi- 
lation of elementary particles via emission of a pair of 
photons. The radiative annihilation leads to a significant 
coupling of a breather to linear waves and brings about 
a variety of resonant and self-oscillation phenomena [16j , 
which can be useful for achieving a coherent superradiant 
emission from coupled systems [161 117] . 



GENERAL RELATIONS 

We consider one-dimensional chains /junctions de- 
scribed by the perturbed sine-Gordon equation: 



Lp — if — aip — sm (y9 — 7, 



(1) 



where <y9 is the phase variables, "primes" and "dots" de- 
note spatial (/?' = dip/dx and temporal (/J = dip/dt deriva- 
tives, a is the viscous damping parameter and 7 is the 
driving (bias) term. In the absence of perturbation terms 
a = 7 = it reduces to the pure SG equation: 

if" ~ (p — soup. (2) 

The soliton in the SG Eq.([2]) is a 2ti phase kink [4]: 

F = 4arctan[exp(a; — ?ii)/v 1 — u^], (3) 



where u is the velocity of the soliton, normalized by the 
speed of light (the Swihart velocity) . The velocity depen- 
dent factor represents the relativistic contraction of the 
soliton when its velocity approaches the speed of light 
u — >■ 1 [1] . This is the consequence of Lorentz invariance 
of the SG equation ([2| . The normalized energy of a static 
soliton M = is Egoi = 8. 



Here the first, the second and the third terms repre- 
sent correspondingly the magnetic/elastic, the Joseph- 
son/potential and the electric/kinetic energies for the 
case of a junction/chain. 

Coupling leads to splitting of the dispersion relation of 
small oscillations into N eigenmodes with different sym- 
metries and propagation velocities [H [H] : 



The coupled sine-Gordon equation 

We assume that a system of N interacting junctions 
can be described by the perturbed CSGE [7j: 



v'l = Aij [(p.j + aLpj + sin Lpj - 7] 



(4) 



1 -25" cos 



iV- 



-1/2 



, (n = l,2,...,iV). (9) 



The slowest mode n — N corresponds to out-of-phase 
(antisymmetric) oscillations in neighbor junctions and 
the fastest, n = 1 to the in-phase (symmetric) oscilla- 
tions in all the junctions [T5] . 



Here i, j = 1, 2, ...N is the junction index and A^j is the 
coupling matrix, off-diagonal elements of which describe 
interaction between different junctions. In what follows 
we will consider the simplest case of nearest neighbor 
interaction, described by a symmetric tridiagonal matrix 
the only nonzero elements of which are: 



Aii — 1, A 



i,i—l 



= A 



i,i+l 



Here S < 0.5 is the coupling strength. The CSGE can 
be also written in the equivalent inverted form 



Ay Vr - 'fij - aifij = simp J - 7. 



(5) 



Apparently, for A^ = 1 it reduces to the perturbed SG Eq. 
(fTl). In case of two coupled chains iV = 2, the system of 
unperturbed CSGE a = 7 = reads: 



1 



S 



;ipj^ - ifl = Sm(y9l 



r¥'2. 



——^ip'l -(p2 = sin (^2 - YT752^i- 



(6) 

(7) 



It is easy to verify by direct application of the Lorentz 
transformation that the CSGE is not Lorentz invariant, 
unlike the SG equation ([2]). 

Physically, the considered type of coupling corresponds 
e.g., to magnetic (inductive) interaction of stacked 
Josephson junctions [7j . In this case space and time in the 
dimensionless equations are normalized by the Joseph- 
son penetration depth, Aj, and the Josephson plasma 
frequency Wp, respectively, the velocity is normalized by 
the Swihart velocity cq = AjWp and 7 by the Josephson 
critical current 7 — I jlc- More details on the normal- 
ization and the formalism can be found in Refs.[9l [TO] . 
As mentioned in the introduction, coupling terms as in 
Eqs. ( |4|5[ ) are also relevant for other objects, like atomic 
chains [2] and magnetic multilayers . 

The energy density of the coupled system is [HI [TO] 



dE{x) 
dx 



= l^'A, 



N 

E 



(l-cos(^,) + 2V''- (8) 



A single soliton in the unperturbed CSGE 

For the solitonic motion with a constant velocity u, 
ip{x,t) — ip{x~ut), the unperturbed CSGE Eq. ([5]) with 
a = 7 = can be written in the simple vector form |10) : 



[A,^ 



u'^E]ip" = simp. 



(10) 



Here E is the unitary matrix. This equation is essentially 
similar to the static CSGE, for which the first integral is 
known 9J. Therefore we can in a similar manner write 
the first integral for the solitonic motion: 



l^',[A7.^ - u^EU ~ S,^i[l - COS(p,;] = C. 



(11) 



For a single soliton the constant C = because at the 
infinity <p>i = p\ = 0. From comparison with the general 
expression for the energy density Eq. ([S]), it is easily 
seen that the soliton energy is twice the magnetic/elastic 
energy 



E,oi = 2£;,„ = / dxp'JA.^p'i, 



(12) 



as is also the case for the soliton in the SG equation [?]. 
The exact soliton solution in the CSGE is not yet 
known. However, an approximate composite soliton solu- 
tion has been proposed, verified by numerical simulations 
and by perturbation correction calculations [H H^ . It is 
represented by a linear superposition of solitonic waves 
Eq. ([3]), corresponding to different eigenmodes. 



^i 



N 



,F„, (* = 1,2,...,7V), (13) 



Fn — A arctan 



exp 



ut 



\n^l-V?/cl 



(14) 



Here i is the junction number, A„ is the characteristic 
length scale of the eigenmode n: 

-1/2 



An — 



1 -25" COS- 



A^ + 1 



(n = l,2,...,iV). (15) 



Note that A„ ^ are eigenvalues of the coupHng matrix Ay 
and coefficients Kn,i are components of the eigenvectors 
of Aij, normalized so that J2n=i'^n,i = 1 in the junc- 
tion containing the soliton and zero in all other junctions 
[TUl FIB] . Thus the soliton consists of a 2tt kink in one 
junction and ripples in all other junctions. The soliton 
shape (coefficients Kn,i) does depend on the junctions 
number and is, for example, different for the soliton in 
the outmost and in the central junctions of the stack. 
Amplitudes of ripples in the neighbor junctions depend 
on the coupling strength and can be significant in the 
strong coupling case S ~ 0.5. The ripples decrease with 
the distance from the soliton both along and across the 
junctions. 

As discussed in Ref. [S] , the static soliton energy in the 
CSGE is larger than that in the single SG equation both 
due to presence of ripples in neighbor chains and recon- 
struction of characteristic length scales Eq.(15). Let's, 
for example, estimate the energy in the simplest case 
N = 2. In this case the multicomponent soliton, Eq. 



(131, becomes 



7V = 2 



2 , A2 = (1 



fl ~ 2 

F-i-F- 
^2 - — ' 



S) 



-1/2 
-1/2 



(16) 



tnat Ji-i 1 — 2 2 
S"^) we obtain 



Substituting those into Eq. ( 12 ) and taking into account 

1/(1 



S'^) and A^l = A:^\ 



5/(1. 



E,oi{N 



2) = 
4[(1 



1-S2 
f 5) 



(17) 
B.92, which is ~ 1.12 
8 for a single junc- 



For S = 0.5 we get Esoi{N = 2) 

times larger than Esoi{N ~ 1) 

tion, in good agreement with numerical simulations [S]. 

The soliton energy increases with TV and saturates at 

~ 3.6Esoi{N — 1) in the strong coupling case S = 0.5 

m- 

The shape of the soliton in the CSGE experiences 



In- 
(|13|), ex- 



strong metamorphoses in the dynamic case [HI HO ■ 
deed, since the soliton components Fn, Eq. 
perience Lorentz contraction at different characteristic 
velocities c„, Eq. ^, the relative shape of the soliton 
does not remain the same as in the static case. When 
u approaches the slowest velocity cat, the correspond- 
ing component n = N gets contracted, while the rest 
of the soliton remains uncontracted. Such the partial 
Lorentz contraction was confirmed by numerical simula- 
tions (HI HH] ■ The soliton survives even at superluminal 
velocity u > cn, however in this case the characteristic 
length Aat becomes imaginary due to the Lorentz factor. 
This corresponds to transformation of the corresponding 
soliton component Fj^ into the out-of-phase plasma wave 
[S] . A similar process of decomposition of soliton compo- 
nents into plasma waves with corresponding symmetries 
occurs when u exceeds any of the characteristic velocities 



Cn [lU]- The phenomenon resembles Cherenkov emission 
from a superluminal particle [9j [TOl [15] with the excep- 
tion that the speed of light is multiple- valued, Eq. (|9|, 
and the dispersion relation is different. 

Multisoliton states in the GSGE are dominated by a 
profound metastability [HlllHISn], i-e-, for given bound- 
ary conditions a large variety of metastable soliton distri- 
butions is possible. Moving solitons interact with linear 
waves, which leads to appearance of geometrical reso- 
nances (standing waves) in finite size systems. Note that 
a soliton in the CSGE can excite all eigenmodes Eq.(|9|, 
which leads to a large variety of geometrical resonances 
El [131. 



Breathers 

Breather is a bound soliton-antisoliton pair. The 
breather solution of the SG equation Eq. (pi) in the 
center-of-mass frame is [3] 



ip — A arctan 



tanz^sin[(cosj^)i] 
cosh[(sini/)a;] 



(18) 



Here < i^ < 7r/2 is determining the breather ampli- 
tude LfBr = 4z^. The breather is oscillating without an- 
nihilation or decay at a frequency lobt — cos u, which is 
always less than the plasma frequency lobt < 1- The so- 
lution Eq. (18) is valid for an infinite system L — cxd. A 



more complicated solution for the finite size system can 
be found in Ref. [2T] . The total energy of the breather 
Esr — 16 sin u is smaller than the energy of two static 



solitons Ebt < 2Es 



16, leading to binding of the 



soliton and the antisoliton. 

Finite dissipation a > leads to decay of the breather 
and facilitates soliton-antisoliton annihilation. The decay 
is primarily caused by the viscous damping of the soliton- 
antisoliton motion. However, minor radiative losses also 
appear [1] . A qualitative change of the wave form takes 
place upon the soliton-antisoliton annihilation. Initially, 
the soliton-antisoliton pair (Av ~ 2tt) shrinks, i.e., the 
maximum separation between the pair Aa; ^ 1/tanj/ ^ 
1 gradually decreases after every collision. At 4i/ < tt, 
the soliton and the antisoliton completely merge and can 
no longer be distinguished. Further decay (reduction of 
v) leads to expansion of the breather. From Eq. ( 18 ) 



it follows that for small i^ <C 1 the size of the breather 
is oc 1/ siniy. Eventually, the breather turns into the 
longitudinal plasma wave with lobt = cos i^ ~ 1 and the 
wave number kx — sin i/ ~ 0. This accomplishes the 
soliton-antisoliton annihilation. 

Breathers play role not only in soliton-antisoliton an- 
nihilation and the opposite (time-reversal) process of cre- 
ation (or penetration) of the soliton [TS] . Breathers also 
interact with traveling waves and external forces [T|. In 
Ref. [ini it was argued that breathers in the CSGE can 
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FIG. 1. (Color online). Snapshots of time evolution of phase distributions upon soliton-antisoliton collision in the unperturbed 
(a = 7 = 0) double-junction CSGE for (a) direct collision of a solton and an antisoliton in the junction-1 and (b) indirect 
collision of a soliton in the junction-1 and an antisoliton in the junction-2. Thick and thin lines represent tpi and (p2, respectively. 
It is seen that both the horizontal (a) and vertical (b) breathers are decaying due to emission of plasma waves, even in the 
absence of perturbations. 



help to pump energy from the external dc-power supply 
(7 term in Eq. Q) into the oscillating travelling waves, 
which leads to appearance of self-oscillations in the dc- 
driven CSGE with 7 > 0. The phenomenon may find 
practical applications for generation of coherent (super- 
radiant) THz sources based on stacked intrinsic Joseph- 
son junctions in high-temperature superconductors jT7] . 

In the CSGE we will consider two distinctly different 
types of breathers (referred to as the horizonal and the 
vertical): 

(i) the "horizontal" breather corresponds to a direct 
collision of a soliton and an antisoliton in the same junc- 
tion. 

(ii) the "vertical" breather corresponds to an indi- 
rect collision of a soliton and an antisoliton in different 
junctions. Unlike the horizontal breather, the vertical 
breather does not annihilate even in the presence of dis- 
sipation, but leads to formation of a stable static soliton- 
antisohton pair [201112]. For TV = 2 the static "vertical" 
soliton-antisoliton pair has an exact antisymmetric solu- 
tion [7^: ipi — —ip2 — F^- The energy of the vertical 
pair 

E^ertiN = 2) = 2 X 8A^i = 16(1 + S)'^/^ (19) 



is smaller than twice the isolated soliton energy Eq.(17l, 
leading to binding of the pair. 



Numerical procedure 

The system of partial differential equations Eq. ^ for 
different N and junction length L is solved numerically 



using an explicit finite difference method (central differ- 
ence in space and time). The spatial mesh size Ax was 
typically 0.025 and the temporal Ai = Ax/10. The ab- 
sence of spurious effects was checked by changing mesh 
sizes and integration times. 

Static solitons and antisolitons Eq. (pi were introduced 
at certain positions at the initial time. The system is then 
given long enough time to relax with a large damping 
factor a = 2. The large viscosity prevents significant 
soliton motion during the transient period. After that 
calculations continued with the desired value of a. The 
time count t = starts from the end of the transient 
period. 

All simulations were made for zero external field 
boundary conditions at x = ±L/2 



difii/dx = 0. 



(20) 



Those boundary conditions are non-radiative, i.e., pre- 
clude energy flow through the edges 123]. In some cases 
dynamic radiative boundary conditions were employed 
(still at zero external field) following Ref. [21]. The 
radiation emission is facilitated by the finite radiation 
impedance Z. For more details see Ref. [24] . 

All presented simulations are done for the strong cou- 
pling case S ~ 0.5, close to the maximum value, rele- 
vant e.g. for atomic scale intrinsic Josephson junctions 
[TU^13| . It was checked that variation of the coupling 
strength does not affect the qualitative presence of the 
effects described below. 
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time approximately at the same time {ti ~ 47). For a 
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FIG. 2. (Color online). Detailed view of the time evolution 
of ip2 (thin lines) for the horizontal breather from Fig. 1 (a). 
The phase ipi (thick lines) is shown close to the moments of 
collisions, marked by arrows. It is seen that emission of two 
wave fronts occurs at every collision: the fast front (marked 
by blue dashed lines) has an in-phase symmetry (pi — <^2 and 
propagates with the fast velocity Ci , the slow front (marked by 
green dotted lines) has an out-of-phase symmetry ipi = ~ip2 
and propagates with the slow velocity cat. 



RESULTS 
Unperturbed soliton-antisoliton dynamics 

Figure 1 shows time sequence of calculated phase pro- 
files ifi (thick lines) and ip2 (thin lines) for the unper- 
turbed (a = 7 = 0) CSGE in a double junction system 
N = 2, Eq. (l7|, for the horizontal (a) and the vertical (b) 
breathers. Initially at i = the soliton and the antisoli- 



single SG equation, the breather Eq. ( 18 1 would continue 



to oscillate without decay with the same periodicity, i.e. 
the subsequent collision would occur at time intervals 2ti. 
This is clearly not the case in the unperturbed CSGE: 

(i) First of all, subsequent collisions occur at smaller 
time intervals. For example, the second collision for both 
breather occurs at ^2 — 59 and t2 — ti ~ 12 much shorter 
than 2ti ~ 94. The third collision for the horizontal 
breather occurs at ty, ~ 68 and ia — ^2 — 9 and so on. 

(ii) Second, the amplitude of the horizonal breather 
decays with time. The soliton and the antisoliton in the 
vertical breather can not change their ±27r amplitudes, 
instead they slow down and eventually form a static pair. 

(iii) Travelling waves are emanating from the breather 
after the collision. 

Figure 2 shows time dependence of (p2 for the horizonal 
breather from Fig. 1 (a) . For simplification the phase (pi 
(thick lines) is shown only at the moments of collisions, 
marked by arrows. It is clearly seen that waves are emit- 
ted from the breather upon the first soliton-antisoliton 
collision ti ~ 47. Two wave fronts can be distinguished: 
the faster (marked by the lower dashed blue lines) prop- 
agate with a constant speed ci = 1.414, and the slower 
(marked by the lower dotted green lines) with the speed 
C2 ~ 0.8165 in agreement with Eq. (|9|. A compari- 
son of (/3i^2 clearly demonstrates that the faster front has 
the in-phase v^i = ^2 and the slower the out-of-phase 
eigenmode symmetry. At the second collision ±2 — 59, 
two new wave fronts are emitted, marked by the upper 
dashed and dotted lines, originated at the second colli- 
sion point {x — Q, t ~ 59). Every subsequent collision 
leads to a similar emission. We emphasize that this hap- 
pens in the absence of dissipation a = 0. Therefore, the 
breather is decaying in the unperturbed CSGE entirely 
due to radiative losses. 

Figure 3 (a) shows time dependence (counted from the 
first colhsion ii) of the total energy Etot (top black line), 
the electric energy E,,, given by the third term in Eq. 
([8]) (middle red line), and magnetic energy Erm given by 
the first term in Eq. (Is]) (bottom blue line) for the case of 
a horizontal breather in the middle i = 5 of A'^ = 10 cou- 
pled junctions with the length L — 100. Calculations are 
made for the unperturbed CSGE a — "f — Q and with- 
out radiation emission at the edges Z = 00. It is seen 
that Etot is conserved because there are no dissipative or 
radiative losses. Maxima in E^ and minima in Em oc- 
cur upon soliton-antisoliton collisions. It is seen that the 
period of collisions is decreasing with time and the mag- 
netic energy Em, related to the breather amplitude, is 
rapidly decreasing after the first collision. This is similar 
to the N — 2 case shown in Fig. 1 (a). At i — ii > i/c„ 
the emitted waves from the breather reflect back from 
the edges and come back to the breather. This leads to a 
very complicated phase pattern consisting of a breather 
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FIG. 3. (Color online). Time dependence of the total Etot, electric Ee and magnetic Em energies for a horizontal breather in 
junction i = 5 of the unperturbed A'^ = 10 CSGE system, (a) Without radiative losses at the edges Z = oo, (b) with radiative 
losses Z = 10^ fl. Panel (c) shows a long-time evolution of -Etot(i) in the presence of radiative losses. A peculiar logarithmic 
time decay is seen. 



and bouncing waves from all A'^ = 10 eigenmodes. 

In order to avoid a complication associated with the 
reflection and bouncing of the emitted waves, we made 
simulations for the same parameters with a finite radia- 
tive impedance Z [21]. In this case the waves partly 
transmit through the edges and leave the system. This 
leads to a decay of propagating waves, except for the de- 
structively interfering out-of-phase mode n = N, which 
can not be emitted (see Ref. [24] for a discussion of emis- 
sion efficiency of different eigenmodes) . 

Fig. 3 (b) presents the results of simulations for 
the same case as in (a) but with the finite radiative 
impedance. It is seen that initially the total energy is 
conserved until t — ti = L/2ci ~ 10 when the fastest in- 
phase (n = 1) wave front reaches the edges. After that 
the energy starts to decay due to radiative losses through 
the edges. At large times the decay of the breather energy 
Etot(t) slows down. Simultaneously, long period beatings 
in Em, due to slow flexural oscillations of the breather be- 
come obvious. 

Do the soliton and the antisoliton completely annihi- 
late upon direct collision in the unperturbed CSGE, or 
do they eventually form a stable non-decaying breather? 
To answer this question we performed long-time calcula- 
tions. In order to avoid possible radiative losses at the 
edges from the tail of the breather itself, we studied even 
longer systems. Fig. 3 (c) shows such simulations for 
L = 300 and the rest of parameters the same as in panel 
(b). A peculiar logarithmic time decay is clearly seen: 



Etot{t)o^Et,t{0)-p\og{t/t,), 



(21) 



where /3 depends on the strength of radiative losses at 
the edges. Thus, unlike in the unperturbed sine-Gordon 
equation, in the unperturbed coupled sine-Gordon equa- 
tion the soliton-antisoliton pair does annihilate upon the 
direct collision even in the absence of dissipation, but the 
annihilation takes an exponentially long time. Therefore, 



for all practical cases the breather would appear stable 
at the time of the experiment, just like the circulating 
current in type-II superconductors [25] . 

The indirect soliton-antisoliton collision, shown in Fig. 
n (b), does not lead to annihilation, but to formation 
of a static bound pair. The vertical breather, produced 
upon the indirect collision decays in a similar way as the 
horizontal breather discussed above. The total radiative 
losses upon the indirect collision is the difference between 
the energy of two isolated solitons and the static vertical 
soliton-antisoliton pair. For the double junction, shown 
in Fig. 1 (b), those are given by Eqs. (17|19): 



AE{N = 2) = 8 (1 - S") 



N-l/2 



(l + 5)-i/2 



(22) 



For S = 0.5, as in Fig. 1 (b), about 27% of the initial 
energy is lost into radiation. 



Soliton-antisoliton dynamics in the perturbed CSGE 

Addition of viscous damping a > perturbation term 
leads to decay of the breather both in the single SG and 
in the CSGE. As seen from Fig. 3, reduction of the am- 
plitude of breather oscillations is accompanied by the 
increment of oscillation frequency. Figure 4 shows the 
amplitude of the breather as a function of breather fre- 
quency for the SG {N = 1) and the horizontal breather 
in the CSGE (A^ = 2 and 10) with different damping a. 
It is seen that in all cases the if Br {oJ Br) dependence fol- 
lows well the theoretical expression Eq. (18 1, shown by 



the dashed line. The scattering of points for the CSGE 
case is due to presence of the strong radiative field, which 
complicates the determination of the breather frequency 
and amplitude, see Fig. 3 (a). 
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FIG. 4. (Color online). Phase amplitude of the horizontal 
breather as a function of the breather frequency for perturbed 
SG {N = 1) and GSGE with iV = 2 and 10 for different values 
of the damping a and for 7 = 0. The dashed line represents 
the breather solution Eq. (fTsl) for the SG equation. 
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So far we considered the case without driving force 7 = 
0. As shown above, in the CSGE soliton and antisoliton 
always annihilate due to either radiative or dissipative 
losses. The driving force 7 7^ replenishes the energy lost 
in the collision and may lead to survival of the solitons 
after the collision. 

In a finite system, a moving soliton will inevitably col- 
lide with the edges. For the zero-field boundary condition 
the collision of the soliton with the edge is equivalent to 
a collision with an image antisoliton 4!. After the col- 
lision the image soliton continues the motion, i.e. the 
soliton is reflected as an antisoliton at the edge [1]. The 
shuttling soliton-antisoliton motion leads to appearance 
of zero-field steps (ZFS) in current-voltage (I-V) charac- 
teristics of Josephson junctions O [26H29] . In this case, 
current is the dc-driving force / = 7 and dc voltage is 
time-average of the velocity. 

Figure 5 shows snapshots of voltage profiles Vi — tpi 
for a single soliton in the junction i = 1 of a double junc- 
tion structure TV = 2 for different driving terms 7 and 
for a = 0.01 and L = 200. The time is counted relative 
to the first collision ti with the left edge. Panel (a) cor- 
responds to a small driving force 7 = 0.01 and a slow, 
subluminal soliton motion u = 0.57 < cn- Before the 
collision, t — ti = —11, the soliton was moving to the left. 
After collision it is reflected as an antisoliton moving to 
the right. Simultaneously, emission of plasma waves from 
both the in-phase and the out-of-phase eigenmodes takes 
place, similar to the unperturbed case in Figs. 1 (a) and 
2. Panels (b) and (c) show snapshots at larger driving 
forces and soliton velocities larger than the out-of-phase 
plasma wave speed cn- Such super luminal soliton mo- 



FIG. 5. (Color online). Instantaneous voltage (velocity) 
V{x) — dipi^2/dt profiles for a driven soliton motion in a 
dissipative a — 0.01 CSGE with N = 2 for increasing driving 
currents (forces) (a) 7 — 0.01, (b) 7 — 0.015 and (c) 7 — 0.03. 
Time is counted with respect to the first collision ii with the 
left edge of the system. It is seen that the fast in-phase and 
the slow out-of-phase waves are emitted upon the collision 
t = ti. In panel (b) the soliton velocity u is close to the 
velocity of the out-of-phase wave, and the corresponding front 
is no longer seen ahead of the soliton. Further increase of u in 
panel (c) leads to profound Cherenkov-type radiation behind 
the soliton. 



tion is accompanied by Cherenkov-type radiation behind 
the soliton (21 [TOl [ISl 130] • A comparison of snapshots at 
t — ti = 144 in panels (a-c) indicates that the in-phase 
radiation front from the collision event is similar for all 
shown soliton velocities, but the out-of-phase front is not 
visible ahead of the soliton when the soliton is moving 
faster than the out-of-phase plasma wave [31] • From this 
it is also clear that it is the edge, rather than the moving 
soliton, that emanates the waves. 



Soliton resonances in finite-size systems 

A shuttling soliton in a finite-size system will periodi- 
cally excite travelling waves at the edges. The emanated 
waves also propagate along the chains and reflect back at 
the edges. The shuttling soliton will interact with bounc- 
ing waves. Resonance will appear if bouncing waves are 
in-phase with the soliton at the edges [35]. In Joseph- 
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FIG. 6. (Color online). Current (driving force) - Voltage (velocity) characteristics of a shuttling single soliton in a moderate 
size systems L — 5 for (a) a single junction (SG, A'' = 1) and (c) a double junction stack (CSGE, N = 2), appearance of a fine 
structure of the zero-field step due to interference with emitted plasma waves is clearly seen. Panel (e) shows the continuation 
oi I — V characteristics at large bias for N = 2 and a — 0.02. Panels (b), (d) and (f) show the corresponding emission powers. 
It is seen that in the double junction system the emission at the velocity matching part of the zero-field step (point A in 
(c,d)) is at minimum, unlike the single junction case (a,b), and the maximum emission occurs at point B, corresponding to 
the in-phase geometrical resonance and the Fiske step in the I — V. Above the ZFS, the system switches to another strongly 
emitting resonance (point C in (e,f)), which represents a breather-type self-oscillation. 
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son junctions this leads to appearance of fine structure 
of Zero Field Steps in I-V characteristics [26l [27] . 

Figure 6 (c) shows calculated dc current-voltage 7 - 
< dipi/dt > (<> indicate averaging in time) character- 
istics for a moderately short L — 5 double junction N = 2 
structure for different damping parameters a. Calcula- 
tions are made for the ZFS mode (1,0), i.e. for a single 
soliton shuttling in the junction i — I. In the second 
junction < d(p2/dt >— 0. The dc ZFS voltage is 

/dipi. 2t:u 

A strong almost vertical velocity-matching soliton step, 
marked as point A in Fig. 6 (c), occurs when the soliton 
velocity approaches the slowest out-of phase eigenmode 
velocity m — > C2. According to Eqs. ^ and (23), for 
S = 0.5, iV = 2 this occurs at < dipi/dt >= 1.026. 
The rapid increase of the driving force 7 at w — > C2 is 
caused both by a partial Lorentz contraction of the F2 
component of the composite soliton [3] and by a rapid 
enhancement of dissipation due to Cherenkov radiation, 
see Fig. 5 (b,c). 

Fine structure is seen at the ZFS below the velocity 
matching step due to resonances between the shuttling 
soliton and bouncing waves emanating upon every colli- 
sion of the soliton with edges. In principle, the soliton can 



interfere and form resonances with any type of periodi- 
cally emitted waves. Those can be waves emitted by the 
fluxon upon passing a defect ^32] or Cherenkov-type emis- 
sion [30] • However, the resonances seen in Fig. 6 (c) are 
different from the previously discussed types. Indeed, we 
consider an ideal system without defects and Cherenkov 
emission does not take place at the corresponding soliton 
velocities, as demonstrated in Fig. 5 (a). 

The observed fine structure of ZFS is due to inelas- 
tic nature of soliton collision in the CSGE, even in the 
absence of perturbations, as shown in Figs. 1 and 2. 
This is specific for the CSGE and is not present in the 
unperturbed SG equation, in which the soliton collision 
is always elastic. Even though, some radiation appears 
in the SG in the presence of dissipation 7 7^ 0, the ef- 
fect is very small [1]. For comparison, in Fig. 6 (a) we 
show ZFS in a single junction case (A^ = 1), calculated 
for the perturbed SG with the same parameters as in 
Fig. 6 (c). The velocity-matching step at m — J' 1 and 
< dif/dt >— 27t/L ~ 1.26 is clearly seen. Unlike the 
CSGE case. Fig. 6 (c), it is entirely due to relativis- 
tic Lorentz contraction of the soliton, Eq. (|3| . Because 
in this case soliton - image antisoliton collision at the 
edges is (almost) elastic, the fine structure is not visible 
(a closer inspection reveals the presence of tiny wiggles 
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at the ZFS). 



Radiation emission 



Coupled systems are interesting from the point of view 
of achieving coherent superradiant emission. In particu- 
lar, THz emission from stacked intrinsic Josephson junc- 
tions in cuprate superconductors at zero applied mag- 
netic field is being actively discussed (THl HZl [231 HI] • 

To estimate the emission from the stack at the ZFS, we 
employed the dynamic radiative boundary conditions, as 
in Ref. [23]. The effective radiative impedance was very 
large so that radiative losses do not affect soliton dynam- 
ics. The emission power from the left edge of the double 
junction stack is shown in Fig. 6 (d). Noticeably, the 
emission at the velocity-matching step A is at minimum, 
despite a large dissipation power P — IV, because at 
point A the soliton is resonating with the (Cherenkov) 
out-of-phase eigenmode n = N = 2, see Fig. 5 (c). Even 
though the oscillation amplitude in each junction is very 
large, destructive interference from the two junctions pre- 
vents emission j24l . This is qualitatively different from 
the single junction A^ = 1 case, shown in Fig. 6 (b), 
in which the maximum emission occurs at the velocity 
matching step and the emission power is correlated with 
the total power IV. 

From Fig. 6 (d) it is seen that the main emission occurs 
at the lower resonance B, which corresponds to the volt- 
age of the in-phase cavity (Fiske) mode (m, n) = (2, 1) 
in the stack pJ^IMj : 



Vrr. 



2L 



(24) 



At this point the shuttling soliton excites the in-phase 
standing wave in both junctions, which leads to construc- 
tive interference and to significant superradiant emission 
outside the stack [24]. Note that the emission power at 
the resonance B is rapidly increasing with decreasing a, 
unlike for the rest of the ZFS. This is a clear indication 
that the geometrical resonance is indeed taking place at 
point J5, because the emission at the Fiske step depends 
on the quality factor of the geometrical resonance and 
increases with decreasing a [24] . 

The appearance of the Fiske step at the ZFS provides a 
clear evidence that the shuttling soliton in the CSGE can 
indeed strongly interact with cavity modes and travelling 
waves due to strong emission upon soliton - image anti- 
soliton collisions at the edges. In a similar manner, the 
soliton also interacts with Josephson oscillations when 
one or several junctions are in the running (McCumber) 
state with < dipi/dt >~ 7/0:. This leads to a large vari- 
ety of resonant states [55] and may lead to self-oscillation 
phenomena at geometrical resonance conditions [TO] . 

Figure 6 (e) shows the continuation of the I-V for the 
same N = 2 double junction with a = 0.02 up to higher 



bias current. It is seen that at 7 > 0.33 the system 
switches from the ZFS (1,0) to another strong resonance 
C before it goes into the Ohmic (free running) state. At 
this resonance both junctions have the same voltage and 
are synchronized in-phase. This leads to a large emission, 
as shown in panel (f). Such resonances were discussed in 
Refs. [TH1[111I13]. They may combine 27r soliton kinks 
with a similarly large amplitude waves, which may be dif- 
ficult to disentangle by just looking at the shapes of phase 
profiles (pi{x). However, we observed that high-order ZFS 
can be clearly distinguished from geometrical resonances 
by comparing the emission frequency: ZFS emit at the 
sub-harmonics of the Josephson frequency [28j or even 
at non- Josephson frequency [29], while self-oscillations 
at geometrical resonances emit at the harmonics of the 
Josephson frequency [16j . 

To conclude, we have studied soliton-antisoliton col- 
lisions in the coupled sine-Gordon equation. It was 
shown that in contrast to the sine-Gordon equation, a 
soliton-antisoliton pair annihilates in the CSGE even 
in the absence of perturbations. The annihilation oc- 
curs via a logarithmic-in-time decay of a breather caused 
by emission of plasma waves. In a dissipative, dc- 
driven case, a similar phenomenon leads to a strong 
coupling between the coupled soliton-antisoliton pairs, 
breathers, and propagating waves, which may lead to 
self-oscillations at the geometrical resonance conditions. 
This phenomenon may be useful for achieving superradi- 
ant emission from coupled oscillators. 
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